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Abstract 

The magnitude of binding energy used in the conventional nuclear theory to ex- 
plain the EMC experimental data, seems to be larger than the one expected. In 
this paper to get sufficient depletion in the binding energy, different oscillator-model 
parameters hoj for different shells and the proton (neutron) structure function that 
have good agrement with experimental data are used. The extracted results for 
4 He, 12 C, 40 Co and 56 i ? e nuclei show that one can get improved results in medium 
x ranges by less binding energy. 
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1 Introduction 



In 1983 the European Muon Collaboration (EMC) reported [1] their measurement of muon 
scattering cross sections per nucleon of iron to deuterium nuclei. The ratios were clearly 
different from unitary and totally unexpected because of high momentum transfer involved 
(Q 2 up to 200 GeV). Several phenomena were proposed to explain theoretically the EMC 
effect and each of them had some successes in limited x (the Bjorken variable) range. The 
conventional nuclear theory proposed by Akulinichev et al, [2] and Dunne et al, [3] seems 
to be able to explain the experimental data in medium and large x ranges. Akulinichev 
et al, used a harmonic-oscillator model with the energy levels taken from the compilation 
of experimental data for the (e, ep) reaction [4]. They showed that the nuclear interac- 
tion of nucleons resulting in their binding and Fermi motion, plays an important role in 
deep-inelastic lepton-nuclear scattering. But the average value of the potential for these 
oscillator levels was about -60 MeV. Dunne and Thomas also used a harmonic oscillator 
model with the levels fitted to measured separation energies [3]. So they failed to get 
sufficient depletion for corresponding single particle energy. If in that time someone had 
used measured separation energies within the naive single-particle model, the discrepancy 
would have gotten worse. The authors of references [2, 3], used only one oscillator-model 
parameter hui for different shells inside the investigated nuclei. So they tried to understand 
the EMC effect in term of the change in the mass scale of a bounded nucleon [5]. Also the 
extracted results by the light front formalism that developed by Miler and Smith [6, 7], is 
not consistent with the binding effect not only in the magnitude of the effect, but also in 
the dependence on the number of nucleons. It should be mentioned that they used not only 
one momentum distribution for any nuclear level but also one nucleon structure function 
parametrization. The used parametrization for nucleon structure function in reference [7] 
considers no difference between protons and neutrons structure function inside the investi- 
gated nuclei. To compare it with other parametrization of nucleon structure functions, it is 
plotted in figure 1 by dash curve. Also the convolution approach represented in reference 
[8], no difference between protons or neutrons momentum distribution on different state 
has been considered. As different shells have different root mean square radius one can 
choose different oscillator-model parameters huj for different shells inside a nucleus. Also it 
should be mentioned that the nucleon structure function used by the authors of references 
[2, 3, 6, 7], is far from the experimental data (see figure 1). So in this paper it is tried to 
extract the EMC results in the conventional nuclear theory for the nuclei 56 Fe, 40 Ca, 12 C, 
4 He, first by considering 'the different oscillator-model parameters hu for different shells 
related to their root mean square radius', and second: 'the free GRV's [9] neutron and 
proton structure functions that have good agreement with experimental data'. 



2 Theoretical formalism 



One can obtain the nuclear structure function of nucleus by main equation [2] 




v=n,p nl 
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where the first sum is over the proton and neutron cases. The second sum is over the 
quantum number of each energy levels. The g% is the occupation number of energy level e n \ 
for proton (N=p) and neutron (N=n). The nucleon distribution inside the nucleus define 
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with z = p n iq/rn N qo and x = Q 2 /2m N qo the bjorken variable for free nucleon. The effects 
of the momentum and energy distribution of the nucleon in the nucleus are included in 
Eq. (2) through <t> n i{p) and e n i, respectively. The magnitude of nuclear binding energy (e n j) 
mainly effects the structure functions in the intermediate x region. The function f A (z) n i 
describes the momentum and energy distribution of nucleons inside nuclei and satisfies the 
normalization rule ^ 

E E / dzgZf N (z) nl = A. (3) 

7V=n,p nl J0 

Akulinichev et al, [2] used harmonic-oscillator nuclear wave function to calculated f A (z) n i. 
For the oscillator-model parameter Hou they used only one parameter for any quantum 
number n, I. We knew in the heavier nuclei the deeper closed shells have different root 
mean square radius and for this purpose, in the harmonic-oscillator model one could use 
[10] (see appendix A) 

<r 2 > ni = ^(2n + / + |), (4) 

where 



hu= (2w + Z + -), (6) 



a 2 = ^f (5) 

and by considering m N = 938.905 MeV, one can find in the natural unit 

41.33 , n , 3 
— = 2n + / + - 

< r 2 > nl 2 
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where < r 2 >J t and Hoj expressed respectively in Fermi and MeV unit. Table 1 shows the 
calculated hu for the nuclei that are investigated here and the table 2 contains the brackets 
that shows the occupation number for different levels. 

The resulting expression for the nucleon distribution f N (z) n i inside the nucleus is [2] 

1 [z)nl 2UJ r(n + / + |)^ ^ tl \t 2 \ U-J 

x r l )r[i + t 1 + i2 + i^(,-i-^]. (7) 

For the F^{x/z) the GRV's LO free proton and neutron structure functions parame- 
terizations used [9] 
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and the X) g extended over all light quarks q — u,d, s. The convolution product and the 
convolution product with [ ] + are defined as usual 
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and 
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a s (Q 2 ) 
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2tt ~ (3 \n(Q 2 /A 2 ) /9 3 [hi(Q 2 /A 2 )] 2 



with p = 11 - 2//3, /3i = 102 - 38//3 and A^ 4 = 0.2 GeV. The LO charm quark 
contribution in Eq. (8) is defined as 
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with a = 1 + 4m 2 /Q and in LO 
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where f3 2 = 1 — (4m 2 /Q 2 )-2(l — A*' = 4m 2 and m c = 1.5 GeV. Figure 1 shows the 

extracted GRV's free proton and neutron structure function that used in Eq. (1). The 
related LO GRV's parton distribution Fortran code can be found in [11]. The calculated 
ratio 
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presented in figure 2. To better express deformation of bonded nucleon structure function 
in comparison with free nucleon structure function, the ratio 

R \ x ) = (17) 

is calculated. Where here N is the number of neutron and the Z is the atomic number 
so A = Z + N. The extracted results for R A (x) are plotted in figure 4. 



3 result, discussion and conclusion 

Figure 1 shows the free nucleon structure function (read curve) that is used by S. V. 
Akulinichev et al, [2]. This is far from the GRV's free nucleon structure function ((F^gjjv + 
^2 grv) /2) that is shown by brown curve [5,7]. It seems, if the corrected nucleon structure 
functions were used in Eq. (1), the calculated Remc( x ) c °uld be improved. It is necessary 
to notice that we can not use the nucleon structure function like read curve even brown 
in figure 1, instead of proton and neutron structure functions for the nuclei like nucleus 
56 Fe because of Z^N. In this paper for extracting all the results, the GRV's proton and 
neutron structure functions are used. The difference between the root mean square of shells 
encourages me to consider different hoo parameters for different shells inside a nucleus (see 
table 1). The plotted results in figure 2, 3, 4 extracted by using: a) Proton and neutron 
structure functions that have good agreement with experimental data, b) The fiuj parameter 
that is related to shell's root mean square for investigated nuclei. By considering the cases 
mentioned above, the binding energies e n \ (see table 2) that are used to obtain the results 
in figure 2, seems to be less than those used in previous published paper. For example, 
compare e n \ = —40 MeV for 56 Fe [2, 3, 16] with e n i = —26 ~ —32 MeV which are used here 
to extract the result for the 56 Fe that is shown in figure 2 (see table 2). The figure 3 shows 
that with this new used Tiuj parameters (see table 1), the contribution of only Fermi motion 
effects specially for 56 Fe (full curve) is similar to the previous work [15] that is plotted by 
dash curves in this figure but we should notice the full curves rise up about 0.05 ~ 0.02 more 
than related Bjorken parameter scale. The reason of this effect is that the hu parameters 
related to the root mean square radius are a bit larger than the fiui parameters for example 
used in 2. We see that the used harmonic oscillator parameters by Akulinichev, et al., are 
related to larger radius than those that the experimental data shows. In the medium and 
large x ranges the neutron structure function is smaller than the proton structure function 
so the extra 4 neutrons in the 56 Fe nucleus cause the EMC ratio of iron comes down about 
3 percent more than the other nuclei even by putting e n i = in medium x range. In the 
figure 4 according to the Eq. (17) the ratio of investigated nuclei's structure function to 
simple sum of equal free protons and neutrons GRV's structure function plotted. One can 
find a bit difference in the extracted results for each nucleus because of different radius in 
the case of neglecting binding energy and considering only Fermi motion effect (dash curve). 
We can trust the results from x ~ 0.2 to x ~ 0.7 [17] But when the effect of binding energy 
added to the Fermi motion effect the difference between the ratio for different nuclei near 
the x ~ 0.7 is decreased (full curves). The magnitude of difference between full curve near 
the x ~ 0.7 is small enough that one interpret this effect as a result of saturation for large 
x ranges. The figure 4 shows the binding effect plays important role in medium x ranges. 
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Appendix A 

To calculate root mean square radius we start from calculated total wave function for 
harmonic oscillator potential 

rj)(r) = R(r)Y(9,<p) 
u(r) 



r 

The equation governing the radial motion is 

% 2 d 2 u(r 



Y(e,<p). (is) 



+ [l(l + l) + V(r)]u(r) = Eu(r). (19) 

in studing bound states, conditions have to be imposed on radial solution u(r): 

lim u(r) — > 

«(0) = 0. (20) 
the normalization of radial wave function leads to the integrals 

l°° R 2 (r)dr = [°° u 2 (r)dr = 1. (21) 

JO JO 

for harmonic oscillator potential i.e., V(r) = \mw 2 r 2 one can obtain the radial laguerre 
equation with the solution [18] 

u nl (r) = N n , t r l+1 e^ r2 L l + 1 / 2 (2ur 2 ), (22) 

where v = mu/2h , N n> i normalization factor and E = hu(2n + I + |). For root mean 
square radius we have 



< r 2 > 
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= [°° K , {r 2 ) l+l e- 2 ^L l + l l 2 {2vr 2 )L l + l l 2 {2vr 2 )r 2 dr. (23) 
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If we put x = 2vr 2 and dr = ^\f^dx then we have 

K, , (^ M ^r ^i e -'L':^(x)L':^(x)dx, (24) 
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where [19] 



So we can write 



J™ x k+1 e- x L k n (x)L k n (x)dx = {n ^, fc)! (2n + k + 1). 
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By using Eq. (25) and putting k — I + | , we have 
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To calculate normalization factor 
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If we put x = 2ur 2 then we have 
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By substituting this normalization factor in Eq. (27) one can find 



< r 2 > nl = (2u) l+1 V8^- 



(n + l + i) 
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= (2n + / + -) 
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Figure captions 



Figure 1: The LO GRV's structure functions for proton (blue) and neutron (pink) at 
Q = 4 GeV without charm quark contribution that is used in Eq. (1). Notice to the 
difference between the GRV's nucleon structure function (brown) and the red curve that 
shows the used nucleon structure function in refrence 2. The experimental data are taken 
from [12, 13]. 

Figure 2: The extracted F^(x) Per nudeonj F\ H [x) Per nuc i eon for A= 4 He, 12 C, 40 Ca and 56 Fe 
from Eq. (1). The F^(x) are taken from GRV proton N=p and neutron N=n struc- 
ture functions according to the Eq. (8) that plotted in figure 1. The used parameters 
{9nii 9nii e m) an d the calculated hu parameter are shown in table 1 and 2. The experimen- 
tal data without error bar are taken from [11, 14]. 

Figure 3: The extracted F^(x) p , I Fl H (x) p . for A= 4 He, 12 C, m Ca and 56 Fe 

O Z V /Per nucleon I A \ / rer nucieon > > 

from Eq. (1) by considering only Fermi motion effect. The full curves were obtained with 
the parameters explained in the caption of figure 2 by putting e n i = 0. The dash curves with 
the same color for the same nucleus were obtained with the parameters hu of Reference [3]. 

Figure 4: The extracted F^/(NFV + NF^) for A=4, 12, 40, and 56 from Eq. (1). The 
used parameters are the same as parameters that are used to extract the results in figure 
1. The dash curves with the same color for the same nucleus shows the extracted results 
by considering only the Fermi motion effect (i.e. e n i = 0). 
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Shell 2 H A He 12 C 28 St A0 Ca m Fe 

Os (2.09,15.35) (1.67,22.23) (1.67,22.23) (1.67,22.23) (1.67,22.23) (1.67,22.23) 

Op (2.44,17.36) (2.44,17.36) (2.44,17.36) (2.44,17.36) 

Od (3.10,15.05) (3.10,15.05) (3.10,15.05) 

Is (3.48,11.95) (3.48,11.95) 

0/ (3.74,13.3) 



Table 1: The brackets contain (< r 2 > nl , TtuS). The oscillator-model parameters huj calcu- 
lated from Eq. (6). 



Shell 2 H 4 He 12 C 28 Si m Ca 56 Fe 

0s (1,1,-1) (2,2,-15) (2,2,-22) (2,2,-20) (2,2,-30) (2,2,-32) 

Op (4,4,-20) (6,6,-20) (6,6,-28) (6,6,-32) 

Od (6,6,-18) (10,10,-26) (10,10,-30) 

Is (2,2,-25) (2,2,-28) 

0/ (6, 10, -26) 



Table 2: The brackets contain (g^, g n ; , e n i(MeV)) for related shell. 
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